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Abstract 

The incorporation of a dividend yield in the classical option pricing model of Black- 
Scholes results in a minor modification of the Black-Scholes formula, since the lognormal 
dynamic of the underlying asset is preserved. However, market makers prefer to work with 
cash dividends with fixed value instead of a dividend yield. Since there is no closed-form 
solution for the price of a European Call in this case, many methods have been proposed 
in the literature to approximate it. Here, we present a new approach. We derive an exact 
analytic formula for the sensitivity to dividends of an European option. We use this result 
to elaborate a proxy which possesses the same Taylor expansion around with respect 
to the dividends as the exact price. The obtained approximation is very fast to compute 
(the same complexity than the usual Black-Scholes formula) and numerical tests show the 
extreme accuracy of the method for all practical cases. 

Key words: Equity options, discrete dividends. 



"The opinions expressed in this article are those of the authors alone, and do not reflect the views of 
Societe Generale, its subsidiaries or affiliates. 

^Arnaud Gocsei is a quantitative analyst in the Equity model validation team at Societe Generale. E-mail: 
arnaud.gocsei@sgcib.com. 

* Fouad Sahel is head of the Equity model validation team at Societe Generale. E-mail: mohammed- 
fouad.sahel@sgcib.com. 



1 



Analysis of the sensitivity to discrete dividends : A new approach for pricing vanillas 



Contents 

1 Introduction 

2 Existing Methods 

3 The method 5 

3.1 Motivations and notations 5 

3.2 Spot/Strike adjustment [6 

4 Numerical tests \8 

4.1 Test on an underlying paying dividends with low frequency 

4.2 Test on an underlying paying dividends with high frequency [ 

5 Conclusion [19 

A Computation of the dividend sensitivities E 

A.l First step: a recursive formula IS 

A. 2 Second step: a martingale argument 2C 

A. 3 Third step: Proof of proposition IA.1I [21 



Gocsei-Sahel 



2 



Analysis of the sensitivity to discrete dividends : A new approach for pricing vanillas 



1 Introduction 

In the classical Black-Scholes framework, we can find in the literature three main ways of 
inserting cash dividends into the model □ : 

1. Escrowed model. Assume that the asset price minus the present value of all div- 
idends to be paid before the maturity of the option follows a Geometric Brownian 
Motion. 

2. Forward model. Assume that the asset price plus the forward value of all dividends 
from past dividend dates to today, follows a Geometric Brownian Motion. 

3. Piecewise lognormal model. Assume that the asset price shows a jump downward 
at each dividend date (equal to the cash dividend payment at that date) and follows 
a Geometric Brownian Motion between those dates. 

Although the first two models lead to a closed-form solution, they are not satisfactory. 
Indeed, the option price obtained in these models is not continuous at dividend dates. 
Moreover, if one considers two options with different maturities T\ < T%, the first two 
models lead to different asset price process dynamics for t < T\, since the dividends paid 
between T\ and Ti are taken into account in one case but not in the other. 

Therefore, it is the piecewise lognormal model which is prefered from a theoritical point 
of view. This paper is dedicated to find a robust pricing proxy for this model. We consider 
an underlying following a Black-Scholes dynamic between dividend detachement dates and 
paying cash dividends at discrete times < T\ < . . . < T n < T, i.e. : 

• for Ti < t < T i+1 : 

dS t = rS t dt + aStdWt, 

• at time Tj : 

Srp-\~ Srp— D<i(^Srp—^ . 

i i i 

where r is the interest rate, assumed constant, W is a standard Brownian motion and 
Di is the dividend policy defined by : 



AOS) 



d if S > d 

s if s < a, 



The cash amounts C\ , . . . , C n are known at the initial date and each d represents the 
dividend cash amount eventually paid at time Tj. 

The dividend policy Di is a liquidator policy as the stock price is absorbed at zero at 
time Tj if < Ci. Consequently, the stock price remains positive. Note that as a practical 
matter, for most applications, the definition of Di(S) when S < Ci has negligible financial 
effects^, as the probability that a stock price drops below a declared dividend at a fixed 
time is typically small. It just ensures the positivity of the price. 

In this paper, we are interested in computing the fair price of the European Call 
Call(<So, K) with strike K and maturity T. Since there is no closed-formula, one should 



x We take here the terminology used in [4]. 

2 This becomes less true when considering large maturities and dividends. 
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recover the price via PDE methods using a finite difference scheme, with boundary con- 
ditions at each T, ensuring the continuity of the price of the Call. This procedure can be 
time-consuming if one considers a maturity T = 20 years and an underlying paying as much 
as one dividend a week. Therefore, when computation speed is at stake, one would prefer 
a fast and accurate proxy for the price. 

We review in the following section three of the existing methods in the literature and 
discuss their limitations. 



2 Existing Methods 

1. Method of moments matching. We approximate the stock price process S by a 
process S with a shifted log-normal dynamic under the risk-neutral pricing measure : 



S t = A + Mexp 



1 '9 'ttt 

--a 2 t + a W t 
2 



The three parameters A, M and a are calibrated so that the first three moments of 
St match the first three moments of St- This method reduces to the pricing of a 
European Call on a modified underlying S, which can be done using the usual Black- 
Scholes formula. This proxy does not work well if the stock pays dividends frequently, 
the maturity is greater than 5 years or the option is deep in-the-money. 

2. In PQ, Bos and Vandermark define a mixture of the Escrowed and Forward models, 
using linear pertubations of first order. They derive a proxy resulting in spot/strike 
adjustment : 

CaR(S ,K) « Call BS 
where Call BS is the usual Black-Scholes function and: 



-Sn 



n 

i=i 

n , 



Ti 



-rTi 



K* 



AT-Ti) 



(1) 

(2) 



i=l 



This proxy works better for at-the-money options and small maturities but results in 
serious mis-pricing for in-and out-of-the- money options and large maturities. 

3. In [2], Bos, Gairat and Shepeleva derive a more accurate proxy than the previous one 
by considering a volatility adjustment : 



Call(5o, K) in Call BS (S* , K,a(S* , K,T)), 



with S* given by ([T]): 




e 2 



-rTi 



2T 



S* ^ 

i=i 



N(a)-N[a-o--^= 



+ 



e 2 

*2 



s 



;c 



N(b) -N[b-2a 



T 

min(Ti,Tj) 

7T~ 
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where N(x) is the normal distribution function and: 

a = ^( log (l) + (, '^ 2/2)r ) ,i,=< ' + ^ 

This proxy will be a good benchmark to test the accuracy of our method presented 
in the following section. 



3 The method 

3.1 Motivations and notations 

Consider Call(£o, K) as a function of the dividends C±, . . . , C n : 

Call(S* ,K) = Call(Ci,...,C re ). 

Although there is no closed- form formula for Call(Ci, . . . ,C n ), we prove in annex IA1 that 
we can still compute explicitely its sensitivities to dividends at the origin. More precisely, 
we have for all k E N and 1 < i\, ...,%< n: 



d fc Call ,a fe Call BS 



(0) = (-1) ^ ^ (flbe-^E^i^, Jf.r) e-^U^-^Ui-W, . 



dC h . . . dc lk 

(3) 

We use this result to derive an accurate approximation of Call(Ci, . . . , C n ). Before 
explaining our method, we need first to introduce some notations. For all functions / of n 
variables x±, . . . ,x n and Va € N, we note T a f the a th order Taylor series at of /: 

an r,t. » 

k=0 ii,...,ife=l K 

We introduce the space A a of functions having the same a th Taylor series at as the 
function Call(Ci, . . . , C n ): 

A a := {f,T a f = T a Call}. 

The order a quantifies how near is f(C%, . . . , C n ) from Call(Ci, . . . , C n ) when the dividends 
are small. This precision increase with a. 

Functions / in A a are naturally good candidates to approximate Call. However, the 
difference Call(Ci, . . . , C n ) — f(C±, . . . , C n ) can be quite big if the dividends are not small 
enough. For instance, for a = 2 and 3, we test the accuracy of the natural choice consisting 
of taking 

/(d,...,C n ) :=T a Call(Ci,...,C n ). 

Figure [37T1 shows the relative error of the price of a European Call when using this approx- 
imation. We assume that the stock pays a fixed dividend C every year and we analyse 
how the relative error varies when we increase C . We can see that both the second and 
third order Taylor series give accurate results when the dividends are small but when the 
dividends increase, they both lead to serious mis-pricing. 
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Relative error (in %) for Call ATM with dividends paid every year 



0.30 
0.20 
0.10 




-0.50 J 

C (dividend amount paid every year) 



Figure 1: Relative Error using the Taylor series for approximation. Numerical parameters: 
S"o = 100, K = 100, r = 3%, a = 30%, T = lOy. 

Therefore, the approximations T a Call, a = 2, 3, are not satisfying. Thus, one need to 
find for a > 2, a function Proxy £ A a , different from T a Call, which gives an accurate 
approximation of Call for all practical values of C\, . . . , C n (not necessarily very small). We 
explain in the next subsection how we determine Proxy E A a . 



3.2 Spot/Strike adjustment 

Like Bos and Vandermark in [1], we search our function Proxy under the form: 

Proxy (d, . . . ,C n ) := Call BS (S* {C x , . . . , C n ), K*(C U . . . , C n )), (4) 

with: 



a n 



S*(Ci, . . . , C n ) =Sp + y~] Oj lv .. i j fe Cj 1 . . . Cj fc , (5) 

k=0 h,...,i k =l 



a n 



K*{C x ,...,C n )=K + Y. E K-,i k Ch---C ik . (6) 

k=0 «i,...,ife=l 

The reason why we perform a spot /strike adjustment is that it allows to recover the exact 
price when the dividends are paid spot or at maturity. 

The coefficients a^,...,^ and ^ lv ..,i fe are calculated recursively. They are entirely deter- 
mined by the two following conditions: 

1- B S^(0) = R g^(0),Vfc<a 1 
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2. We impose our proxy to satisfy the Call-Put paritjH: 

n 

Call BS (S*,iT) - Put BS (S*,iT) =S - Ke~ rT -^de- rT \ (7) 

i=i 

n 

i.e. S* - K*e~ rT =S - Ke~ rT - ^ Qe^. (8) 

i=i 

Let's detail the calculus: 

• Computation of a* and the equality 8 gg xy (0) = ^P(O) reads: 

7V(di)ai - e- rT N(d 2 )bi = -e~ rT > N(d(Ti)), (9) 

where : 

di =^ 7 =(ln(So/^) + (r + a 2 /2)T), 

d 2 =di - <t\/T, 
d(t) =d x - S=t, 0<t<T. 



The differentiation of ([8]) writes: 

Oi - e~ rT bi = -e~ rT \ (10) 
Solving the linear system ([9^- (fTU]) gives: 

_ rTi N{d{T l ))-N{d 2 ) 
N(d 1 )-N(d 2 ) ' 

= r( T- T< ) ^(rfi)-^(rf(Ti)) 
N(d 1 )-N(d 2 ) ' 

Computation of etjj and bij: the equality g^g^ 7 (0) = dC G dC W anc ^ ^ wo succes i ve 
differentiations in ([8]) give the following linear system: 

flWoij - e- rT N(d 2 )b hj =0, (11) 
Oij - e-^feij =0, (12) 

where: 

_ 9 2 Call, n , 5 2 Call BS . _ . , , . <9 2 Call BS , „ „. ^<9 2 Call BS , „, 

P = dc . dCj (°) ~ Q * a i g 5 2 ( g o,-*0 " OA' + a j ft i) g g g g fti&j gg2 (5 ,-g). 

After some direct computations, we obtain: 



7 

A — rT 



a + 6 jV(d(Ti)) + N(d(Tj)) + cN(d(Ti))N(d(Tj)) + cfe CT J *JV (d(T* + I))) 



3 The right term in equation Q is not rigorously exact since e rT _E[S , T] is not equal to So — Y17=i rT ' 
but the two quantities are very close. 
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with: 



ry =aSVTN'(d 1 )(N(d 1 ) - N(d 2 ) 



N(d 2 )N (di) - N(dt)N (d 2 



N (dx) - N (d 2 )J [N{d 2 )N (dx) - N{d x )N (d 2 ) 

jv(di)-jv(d2;" 



d=N (di)[N(d 1 )-N(d 2 ) 



• Computation of . i fc and 6^ ... j, , fc > 3: the previous method can be reproduced 
recursively. Knowing all the values Oji,...,j m and bj lt ,,,j m , m < k — 1, we obtain a^...^ 
and fej lj ... i i fe by solving a linear system of the form : 



^Cail 13 ^ _ a fe Call BS . 

g 5 fc — (^Oj-KJOi!,...,^ + — Qx^{So,K)b h) ... tik =u, 

a h,...,i k ~ e &u,...,j fe =0. 

We have presented a simple and general method to derive a function Proxy in A a for 
any a £ N. As for the order a that we choose effectively for our tests, the second order 
computation is a good choice for performance and accuracy. Before presenting the numerical 
results, we recall some desirable properties of our second order proxy 

1. fast computation, even when one considers a large number n of dividends. 

2. recovery of exact price when all dividends are paid spot or at maturity. 

3. arbitrage free with the Call-Put parity. 

4. guarantee of the continuity of the Call price at dividend detachement dates. 

5. accuracy for all practical configurations, even for the extreme cases (deep in-the- 
money-option with large maturity and high frequency of dividends) for which the 
already existing methods of the financial literature might lead to serious mis-pricing. 



4 Numerical tests 

4.1 Test on an underlying paying dividends with low frequency 

We test the accuracy of our proxy on a stock with the following parameters: Sq = 100, 
r = 3%, a = 30%. We suppose that the stock pays a dividend of 3 in the middle of 
every year. We compute the Call price with strike K £ {50, 75, 100, 125, 150, 175, 200} and 
maturity T G {5, 10, 15, 20} using four methods: 

1. the finite difference method, 
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2. the method of moments matching. 

3. the spot/vol adjustment of Bos, Gairat and Shepeleva[2], 

4. our proxy with spot/strike adjustment given by ©-([S]), 

Remember that no approximation is made in the finite difference method. The results 
are given in the following tables. 
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Maturity=5 years 



K/S 


0.5 0.75 1 1.25 1.50 1.75 2 


Price: 

FD (exact price) 
Method of moments 
Proxy BGS 
Proxy GS 


47.14 33.85 24.42 17.79 13.12 9.79 7.39 
47.17 33.87 24.42 17.78 13.10 9.77 7.38 
47.11 33.84 24.42 17.80 13.13 9.81 7.41 
47.14 33.85 24.42 17.79 13.12 9.79 7.39 


Relative error (in%): 

Method of moments 
Proxy BGS 
Proxy GS 


0.07 0.06 0.01 -0.05 -0.11 -0.16 -0.20 
-0.05 -0.05 -0.02 0.03 0.10 0.17 0.24 
0.00 0.00 0.00 0.00 0.00 0.00 0.00 
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Maturity=10 years 



K/S 


0.5 0.75 1 1.25 1.50 1.75 2 


Price: 

FD (exact price) 
Method of moments 
Proxy BGS 
Proxy GS 


46.85 38.21 31.66 26.58 22.56 19.34 16.71 
47.07 38.38 31.77 26.64 22.59 19.34 16.69 
46.65 38.08 31.59 26.57 22.60 19.41 16.81 
46.85 38.21 31.66 26.58 22.56 19.34 16.71 


Relative error (in%): 

Method of moments 
Proxy BGS 
Proxy GS 


0.49 0.45 0.35 0.23 0.11 -0.01 -0.13 
-0.43 -0.35 -0.21 -0.04 0.15 0.35 0.55 
0.00 -0.01 -0.01 -0.01 -0.01 -0.01 -0.02 
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Maturity=15 years 



K/S 


0.5 0.75 1 1.25 1.50 1.75 2 


Price: 

FD (exact price) 
Method of moments 
Proxy BGS 
Proxy GS 


46.47 40.48 35.73 31.85 28.63 25.91 23.59 
47.12 41.04 36.18 32.21 28.91 26.13 23.75 
45.79 40.01 35.43 31.70 28.60 25.98 23.74 
46.49 40.49 35.73 31.85 28.63 25.91 23.59 


Relative error (in%): 

Method of moments 
Proxy BGS 
Proxy GS 


1.41 1.38 1.27 1.13 0.98 0.82 0.66 
-1.47 -1.19 -0.85 -0.49 -0.12 0.24 0.60 
0.02 -0.01 -0.02 -0.03 -0.03 -0.04 -0.04 
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Maturity=20 years 



K/S 


0.5 0.75 1 1.25 1.50 1.75 2 


Price: 

FD (exact price) 
Method of moments 
Proxy BGS 
Proxy GS 


46.02 41.74 38.22 35.26 32.72 30.51 28.57 
47.35 42.95 39.30 36.21 33.55 31.24 29.20 
44.33 40.47 37.30 34.63 32.33 30.32 28.55 
46.10 41.76 38.23 35.26 32.71 30.50 28.56 


Relative error (in%): 

Method of moments 
Proxy BGS 
Proxy GS 


2.89 2.90 2.82 2.69 2.54 2.38 2.21 
-3.71 -3.07 -2.44 -1.82 -1.23 -0.65 -0.09 
0.14 0.03 -0.02 -0.04 -0.05 -0.06 -0.07 
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4.2 Test on an underlying paying dividends with high frequency 

We now check the accuracy of our proxy on an underlying paying dividends every week. 
This situation occurs when considering an index like S&P 500 or Eurostox 50. We take 
the following parameters: So = 3000, r = 3%, a = 30%. We suppose that the stock pays 
a dividend of 2 every week. We compute the Call price with strike K such as K/Sq G 
{0.5, 0.75, 1, 1.25, 1.5, 1.75, 2} and maturity T E {5, 10, 15, 20}. 
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Maturity=5 years 



K/S 


0.5 


0.75 


1 


1.25 


1.50 


1.75 


2 


Price: 
















FD (exact price) 


1359.87 


972.67 


699.65 


508.71 


374.45 


279.07 


210.47 


Method of moments 


1361.05 


973.29 


699.70 


508.39 


373.97 


278.54 


209.98 


Proxy BGS 


1358.88 


972.02 


699.52 


508.99 


375.00 


279.75 


211.21 


Proxy GS 


1359.87 


972.69 


699.68 


508.73 


374.47 


279.07 


210.47 


Relative error (in%): 
















Method of moments 


0.09 


0.06 


0.01 


-0.06 


-0.13 


-0.19 


-0.23 


Proxy BGS 


-0.07 


-0.07 


-0.02 


0.05 


0.15 


0.25 


0.35 


Proxy GS 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 
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Maturity =10 years 



K/S 


0.5 


0.75 


1 


1.25 


1.50 


1.75 


2 


Price: 
















FD (exact price) 


1319.62 


1075.07 


890.03 


746.82 


633.81 


543.14 


469.40 


Method of moments 


1327.05 


1080.52 


893.49 


748.67 


634.44 


542.91 


468.55 


Proxy BGS 


1311.37 


1069.84 


887.68 


746.80 


635.57 


546.23 


473.42 


Proxy GS 


1319.68 


1075.04 


889.96 


746.72 


633.69 


543.02 


469.25 


Relative error (in%): 
















Method of moments 


0.56 


0.51 


0.39 


0.25 


0.10 


-0.04 


-0.18 


Proxy BGS 


-0.63 


-0.49 


-0.26 


0.00 


0.28 


0.57 


0.86 


Proxy GS 


0.00 


0.00 


-0.01 


-0.01 


-0.02 


-0.02 


-0.03 
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Maturity =15 years 



K/S 


0.5 


0.75 


1 


1.25 


1.50 


1.75 


2 


Price: 
















FD (exact price) 


1287.50 


1122.13 


990.74 


883.63 


794.68 


719.49 


655.20 


Method of moments 


1308.65 


1140.08 


1005.32 


895.16 


803.52 


726.17 


660.09 


Proxy BGS 


1258.42 


1102.31 


978.74 


877.99 


794.07 


723.01 


662.03 


Proxy GS 


1288.47 


1122.33 


990.66 


883.42 


794.31 


719.10 


654.79 


Relative error (in%): 
















Method of moments 


1.64 


1.60 


1.47 


1.31 


1.11 


0.93 


0.75 


Proxy BGS 


-2.26 


-1.77 


-1.21 


-0.64 


-0.08 


0.49 


1.04 


Proxy GS 


0.08 


0.02 


-0.01 


-0.02 


-0.05 


-0.06 


-0.06 
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Maturity=20 years 



K/S 


0.5 


0.75 


1 


1.25 


1.50 


1.75 


2 


Price: 
















FD (exact price) 


1260.33 


1144.53 


1049.11 


968.59 


899.43 


839.22 


786.23 


Method of moments 


1303.36 


1183.64 


1083.92 


999.27 


926.36 


862.79 


806.82 


Proxy BGS 


1184.43 


1088.21 


1008.55 


940.91 


882.41 


831.09 


785.57 


Proxy GS 


1264.53 


1145.94 


1049.44 


968.43 


899.04 


838.71 


785.66 


Relative error (in%): 
















Method of moments 


3.41 


3.42 


3.32 


3.17 


2.99 


2.81 


2.62 


Proxy BGS 


-6.02 


-4.92 


-3.87 


-2.86 


-1.89 


-0.97 


-0.08 


Proxy GS 


0.33 


0.12 


0.03 


-0.02 


-0.04 


-0.06 


-0.07 
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5 Conclusion 

We have presented a new approach to deal with cash dividends in equity option pricing 
in a piecewise lognormal model for the underlying. Our method relies on the derivation 
of an analytic formula for the sensitivity to dividends of a European option. We obtain 
a closed-form formula for a European Call which gives both very accurate results for all 
practical cases. 

A Computation of the dividend sensitivities 

Consider a European option of maturity T with payoff h(Sx), with S the stock price follow- 
ing the piecewise lognormal dynamic presented in the introduction. We note its fair price 
at time 

n(5o, T, Ci, . . . , C n ) 

We denote: 

n B5 (s ,r) 

the fair price of the option if S does not pay dividends. The partial derivatives 

dc il ...dc ik {So ' T '°'-' 0) 

are related to the usual Black-Scholes greeks by the following formula: 
Proposition A.l For k G N and 1 < i\ < . . . < < n, we have: 

(So,T,0, . . . ,0) = (-1)*^ (S e-^U^, T ) e-r^T^^-m. 



dC h ...dC ik K u ' ' ' v ' dS k 

Follows a proof of this formula. 

A.l First step: a recursive formula 

We introduce some notations: 

• We define the natural filtration (J r t)t>o associated with the brownian motion W. We 
suppose the filtration right continuous. 

• We define for all < h < t 2 : 

X, t — P (r-° 2 /m2-ti)+v(Wt 2 -W tl ) 

• We denote <f)(So,S,t) the log- normal density associated with the variable SoXq^ 

• We define the functions of n + 1 variables (/i«)o<i<n such as: 

hi{S Ti ,C u . . .,C n ) := e-^ T - T ^E[h(S T )\T Ti }. 
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For the sake of simplicity, when there is no confusion, we will simply denote hi{S) instead 
of hi(S, C\, . . . , C n ). Note that we have n(So, T, C±, . . . , C n ) = ho (So). We can compute 
the functions hi recursively beginning with h n : 

h n (S) = U BS (S,T-T n ), 

and by conditioning, \/i < n — 1: 

h t (S) =e- r ^- T ^E[h l+1 ((SX t ^ u+1 -C l+1 ) + )\T Ti ] 

=e -r(T 1+1 - Tl) h i+1 (S i+1 -C i+1 )<f>(S,S i+1 ,T i+1 -T i )dS i+1 , (13) 

Now, we show how these relations allow us to compute recursively the partial derivatives: 

d k U 



dC n . . . dc ik 



(5 ,0,...,0), 



for 1 < i < n, I; £ N* and 1 < i\ < . . . < i^ < n. First, note that a direct application of 
the theorem of differentiation under the integral sign in the relation (|13|) proves that the 
functions hi, < % < n are infinitely differentiable. Then, using the markov property of the 
log-normal densities: 



4>(Si, Si+i, ti)4>(Si+i, Si+2> U+i)dSi+i — 4>(Si, Si+2,U + U+i), 







we obtain: 

d k h h 



(So,0,...,0)=-e- rT nE 



dC h . . . 8C lk 



dsoc l2 . 



Qg-( s oXo^T ix ,0, ■ ■ ■ ,0) 



= -e~ rT n / * (S il ,0,-~,0)<KS o ,S il ,T il )dS il . 

J 12 ' ' ' f'k 

(14) 

This relation will be very useful for a recursive proof of proposition I A. II since it reduces by 
one the number of differentiations with respect to the dividends. 

A. 2 Second step: a martingale argument 

The proof of proposition I A. II relies heavily on this simple but crucial lemma: 

Lemma A. 2 Consider a process following a Black-Scholes dynamic St = Soe^ r ~ u2 l 2 ^ t+aWt , 
< t < T. Then, for all integer k > and for all real number a > 0, the process: 



Z, :-- 



Qk U BS 



H ' 8S k 
is a martingale. 

The following corollary is easy to derive. 



— ( S t e ka2{t ~ a \T - t) e (k-^)(r+ka 2 /2)t 
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Corollary A. 3 For k£N*,0<t<T and a > 0, we have: 



E 



td k Tl BS 

dS k 



S t e k<j2{t - a \T -t 



k-nBS 



d k U 



dS k 



S e- kr72a ,T)e 



-(k-l)(r+ka 2 /2)t 



(15) 



Proof of lemma \A.£k The drift of Z% is: 



D (k-l){r+ka 2 /2)t 



(k - l)(r + ka 2 /2) 



Qkl\BS gk+ljjBS 

dS k dtdS k 



+ (r + ka*)S t e^ d -^ + ^S^^ 9 -^ 



BS 

ds k +~ 



(16) 



where all the derivatives in the last formula are evaluated in (Ste^ 2 ^ a \T — t). Remember 
that U BS satisfies the Black-Scholes PDE: 

dU BS 1 2c2 d 2 U BS 8U BS ^ n 
dt~ + 2 ~dS^ + rS ^T ~ rU = °" 

Now, differentiate k times this equation with respect to S: 



gk+i^BS d k Ii BS 
dWS k '2- BS** +(r + ^- g5 ^- + (*-l)(r+i^/2)- ggs - = 0.(17) 



d k+l U BS 1 d k+2 U BS 

■ + -a S 2 



One immediately checks that the term in bracket in (|16j) is equal to the left term in ([T 
evaluated in (Ste ka (* _a ),T — t), and therefore is equal to 0. 



□ 



A. 3 Third step: Proof of proposition IA.1I 

We argue by recurrence on the number k of differentiations with respect to the dividends: 



If k = 0, the proposition is trivially true as it simply says: 

n(s ,T,o,...,o) = n B5 (s,T). 

Now, suppose that the property is true at rank k. We want to prove that it remains 
true at rank k + 1. We have by relation (fbT 



d k+i u 



-(S ,0,...,0)=-e- rT 'iE 



dC h ...dC ik+1 
By hypothesis of recurrence, we have: 



d k+1 hi 

dSdc i2 ...dc ik J SoXo ^ ' °' • • • ' 0) 



(18) 
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8 k h- 

n (5,0,. ..,0) 



dC- l2 . . . dCi k+1 

= Jse-^E^^-r.Jy.y.j x e -rE^ a 1 (r 4 ,-T 4l )-<7 a E^i(9-2)(Ti 9 -T il ) j 

= ( 1 ) fc ~^F~ [Se {k+1)(T2{T ^- a \T - Ti^j x e [( fc + 1 ) r +l fc ( fc - 1 ) CT2 ] T ! i-( fc + 1 ) ra - CT2 Ej; 3 1 (<?-2)T i9 
where we set: 

fc+i 



a = 



We differentiate with respect to 5: 



<9S3C 42 . . . 8C ik+1 



(5,0,..., 0) 



dS k 

Inserting this formula into (|18p and using corollary IA.31 we get: 
d k+1 U 

= (- 1 )" +l5 ^R^ (Se-< M > A ,T) exp ( -(* + l)(r + «r> - «r 2 ][> " 2 > T <<. 

\ 9=3 



□ 
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